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Two general problems arise when interpreting the recent cosmic ray data as signals of Dark Matter 
(DM) annihilation: (i) the required cross section is too large by O(IOO), and (ii) the annihilation 
products seem to be mostly leptonic. We propose to address these two problems by assuming that 
the couplings of DM to leptons grow with time. This can be achieved by a dynamic localization 
of DM in extra dimensions. A possible outcome of this proposal is a time (red-shift) dependent 
annihilation signal, in terms of strength and dominant final states. 



0\ 

o 
o 



Oh' 
I ■ 

Oh' 

CD 



(N 
> 

o 
m 

o 
o 



X 



The evidence for the presence of Dark Matter (DM), 
a weakly interacting cosmic ingredient which is stable 
on Hubble time scales, has been gaining strength over 
the past several decades. A confluence of data from rota- 
tion curves of galaxies, observations of Cosmic Microwave 
Background (CMB), and other large scale surveys of the 
universe has established that DM makes up about 25% 
of the energy budget of the cosmos [|] . No plausible can- 
didate exists within the Standard Model (SM) of particle 
physics for DM, making it one of the most compelling 
pieces of evidence that require an extension of our cur- 
rent microscopic understanding of Nature. Among the 
possible candidates for DM, Weakly Interacting Massive 
Particles (WIMP's) are generic in models that address 
the stability of the Higgs potential against large quan- 
tum correction, or the "hierarchy" problem. Here, we 
will focus on this class of particles. It turns out that the 
weak scale of 0(100) GeV sets the correct level of DM 
annihilation cross section in the early universe, resulting 
in an acceptable thermal relic abundance. 

The search for DM is carried out in several ways, 
encompassing direct laboratory detection and indirect 
methods which rely on the annihilation signals of DM. 
Of these latter signals, possible hints in cosmic ray data 
have recently received a great deal of attention. In par- 
ticular, the electron flux anomaly observed by the ATIC 
experiment Q and the excess positron fraction measured 
by the PAMELA satellite |^ have been subject of a flurry 
of activity. There are other experiments that have pre- 
sented signals, potentially consistent with observations of 
DM. However, we will only concentrate on the aforemen- 
tioned cosmic ray data in the 0(100) GeV range and will 
not attempt to account for any other signals. In partic- 
ular, we will focus on two apparent problems that arise 
in interpreting the measurements as pointing to WIMP 
annihilation (we note that purely astrophysical e^e~ pro- 
duction by nearby pulsars has also been considered as an 
explanation of these data in, for example, Refs. 0, ||; 
though this remains an interesting possibility, we will 
discuss only a WIMP-based solution here): (i) the re- 
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quired annihilation cross section is too large by O(IOO) 
1^, if one assumes that DM density was set in the early 
universe, when similar annihilation processes froze out, 
and (ii) assuming that DM is a weak scale particle, it is 
puzzling that the cosmic ray data do not reveal an excess 
in the anti-proton flux. 

Recently, (i) and (ii) have been addressed, for example, 
by postulating a light degree of freedom that couples to 
DM and mediates a long range force j7j , or in the context 
of the Stueckelberg mechanism Q . These are interesting 
possibilities that deserve attention, and they may also 
lead to distinct signals at the Large Hadron Collider 
^. However, we note that, in these proposals, the mass 
scales of the WIMP and the mediator have to be chosen 
judiciously such that the desired effect is obtained. The 
resolution to question (ii) above has also been discussed 
in models of leptophilic WIMP 0. 

In this work, we propose that both (i) and (ii) can be 
addressed by assuming a time-varying leptophilic cou- 
pling of the WIMP. In our proposal, we do not constrain 
the early universe couplings of DM to be strictly lep- 
tophilic and only require this to be a late-time feature. 
In fact, in our framework, the primordial couplings of 
DM to the SM can even be baryophilic and leptophobic. 
The main requirement is that DM couples predominantly 
to or today and that these or other SM couplings 
were somewhat smaller at freeze out. This is the ba- 
sic idea that we present here. However, we will discuss 
a specific realization of this scenario in the context of 
models with extra dimensions [pd]| . One may place the 
SM field content in some of the extra dimensions. In 
particular, fermions can be localized at different points 
along the compact dimensions, giving rise to suppressed 
coupling among various fermions, based on their location 
|l2| . Here, we note that DM with time varying proper- 
ties has been considered in Rcf. ||l^, and more recently 
in Ref. [ |T^ . Dynamic fermion localization, in the context 
of baryogenensis, has been studied in Refs. [|l5[ |l^, [iTt . 

On general grounds, any time variation in the DM cou- 
plings can lead to red-shift dependent annihilation sig- 
nals. In our proposal, one would expect that lepton final 
states become more dominant at smaller red-shifts. Of 
course, this depends on the exact rate of time variation in 
the leptophilic couplings. Depending on the specific de- 



tails, one could end up with the early halo annihilations 
that are dominated by hadronic final states; here, galax- 
ies at high red-shift may give distinctly different signals. 
In certain cases, we could also have a situation in which 
dark matter annihilation becomes inefficient over a non- 
negligible cosmological epoch after the early freeze-out. 
Finally, in the extra dimensional scenario that we con- 
sider, the hierarchy problem motivates compactification 
radii of 0{\/TeV) or larger, with the usual attendant 
collider signatures. 

We will adopt the model presented in Ref.||l^, as 
a concrete example, and for simplicity, will focus on 
the fermion content in one extra dimension y of size 
R > 1/TeV. The field content considered here will suffice 
to demonstrate our main ideas. We will make no attempt 
at making the setup fully realistic, lest the key features 
get lost in unnecessary details. Further components can 
be added to our model in a more expansive work. 

Let X denote a singlet fermion and S a singlet scalar. 
Here, x g^ts localized along the extra dimension by its 
coupling with the bulk profile of S, via the 5D Lagrangian 



^5 = J d^xdyxii^ + (y) 



(1) 



where ^ = 7'' 9^ +7^95 is the 5D derivative and ms^ is a 
bulk Dirac mass for x- Following Ref . |^ , we assume the 
simple harmonic oscillator profile S^{y) — How- 
ever, this profile ultimately reaches constant values, giv- 
ing rise to a wall of thickness A < R [|2| . The resulting 
fermion wave-function for the zero mode of x is a 
Gaussian centered at y^ — 171^^/(2^'^) 
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Let be a 5D scalar and a doublet of the SM SU (2) 
that does not develop a vacuum expectation value (vev); 
{(j)) = 0. We will assume that both x and (j) are charged 
under a Z2 symmetry. The following 5D operator yields 
the leading interaction between x and (p 



Oi = A0%/Ax0i, 



(3) 



where A^vA is a 5D Yukawa coupling [g_8[; L is a lepton 
doublet of the SM. One can easily localize the fermions of 
the SM, using straightforward extensions of Eq. (|l|). We 
note that in order to identify % as a DM candidate, we 
will assume that the zero mode masses satisfy > m^. 
We may include a term m'^(j)^(j) in the 5D action to endow 
(j) with mass > 10^ GeV. The mass of the zero mode 
for X may be induced in a variety of ways. For now, we 
assume that a Majorana mass ra^x^Cx, with C the 5D 
charge conjugation operator, can be induced in the bulk. 
Given the typical expectation rriy ^ 100 GeV, this is a 



small mass compared to other scales and we treat it as a 
perturbation. 

We will further assume that the profile of the SM Higgs 
doublet is flat over the width of the wall. Bulk masses 
TO5/ of the SM fermions allow one to localize them at 
various points yf along the fifth dimension. Given the 
exponential nature of the fermion wave-functions, a wide 
range of 4D masses m f can be obtained without recourse 
to small Yukawa couplings 

Here, vew ^ 10^ GeV is the Higgs vev and s is the 
separation of the SM SU{2) doublet and singlet fields 
along the fifth dimension; we have assumed 0(1) Yukawa 
couplings and a separate localizing profile for the SM 
fields, characterized by /i/, which is not strictly necessary. 
In this way, the lightest fermions are typically the farthest 
separated. We will mainly focus on lepton localization. 

The proximity of x to the first generation lepton dou- 
blet Le, at the present epoch, ensures a significant cou- 
pling only to these fermions. If the the heavier charged 
lepton doublets i/j,r are mildly separated from Lg, they 
would not couple to the DM candidate x significantly. 
This feature can naturally account for xx ~* e^e~ as 
a preferred (^-mediated present-day-annihilation channel. 
Note that in this particular case, the dominant coupling 
is to the left handed fields (e, i^e)L, and hence annihilation 
into Ve pairs (or other neutrino flavors, corresponding to 
different dominant charged lepton final states) can be 
comparable to that in the e'^e~ channel. This could lead 
to possible experimental signatures in astrophysical high 
energy neutrino experiments; see, for example, Ref. 
for a discussion. 

With a roughly flat zero mode profile for cf), the an- 
nihilation of X is dominantly governed by the coupling 
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where uj^ = / + /i^) and s^e = \yx - Ve\ is 

the separation between the Gaussian centers; different 
localizing scalars for Lf. and x are assumed. 

If we keep the masses of the particles fixed, the annihi- 
lation cross section of x is of order (ctow) ~ ^^e/'^^' ^i^^^ 
™x ^ ™0- Hence, getting a factor of 0(100) enhance- 
ment in the cross section at late times roughly requires 
9xe ~* 3(7^e, approximately one e-folding, which only 
amounts to an 0(1) shift of the exponent in Eq. (^). Be- 
fore discussing the required shift, let us simplify Eq. (^) 
by assuming that all fermions are localized by the same 
profile and /i^ = /^e = We then have, 
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s^g — 2/^^. Since we are only 



Thus, we need s^g 

interested in the relative distance between the fermions, 
let us assume, for purposes of demonstration, that y^ = 0. 
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Then, we find that the new position of x should be at 
y'x = "^5x/(2m^)' where 

Note that accommodating a realistic set of masses in the 
SM typically requires /is to be of order a few . Hence, 
one could typically expect m^f / ^ ^ 1, leading to m'^^ « 
777-5^(1 — 'ijj? / . Thus, we estimate 

l^»T^5xl^M- (7) 

This relation quantifies the shift in the value of ms^, in 
terms of the localizing potential parameter fi, required to 
account for both the primordial and present day annihila- 
tion cross section of x- We have implicitly assumed that 
both rates are dominated by the interaction in Eq. (||). 
This in turn implies that other lepton doublets did not 
have a major overlap with x, during the evolution of g^g. 

Let us now estimate what Eq. (|^) entails in some moti- 
vated scenarios. Given that A would set the typical mass 
scale of Kaluza-Klein (KK) modes of SM particle trapped 
within the wall, we expect A"^ > 10^ GeV Q; we will 
take A^i - 1 TeV. Following Ref. we may then 

expect ^l ~ lOA^i ~ 10 TeV. So far, we have assumed 
the 5D mass parameter for fermions, and in particular 
X, to be a Lagrangian parameter. However, in general, 
we may postulate that this mass is set by the vev of a 
5D scalar $ (note that this field is not that in Ref. ||T^ 
responsible for localizing the fermions). Assuming a con- 
stant rate of change of ($) over the cosmological time 
scale to ^ Hq^, where Ho ~ lO^'^'^ eV is the present day 
Hubble constant, and (5$ ^ /i we find 

$2 ^ {^Ho f - 10-4° eV^ (8) 

The inferred cosmic dark energy density is 
0(10"") eV* 0. So, we do not expect this rate 
of change for $ to have a significant effect on the 
cosmological evolution of the universe. Nonetheless, 
obtaining such a tiny rate of change for $ is fraught 
with difficulties, as this field must have an extremely 
small mass, m$ < Hq, so as to vary over cosmic time 
scales (a familiar feature in quintessence models). As 
we need $-x coupling of C(l), maintaining such a 
fine-tuned potential poses a theoretical challenge. Also, 
a nearly massless scalar with significant couplings to 
DM could cause violations of equivalence principle, 
which is potentially another undesirable feature. Here, 
we instead relegate the variation of to a temporal 
change in the coupling ^ of <i> to x and hold ($) fixed. 
We will consider cases where ^ depends on dynamical 
moduli in the gravitational sector. The couplings of 
such moduli are typically suppressed by the Planck 
mass Mp ^ 10^^ GeV, which in turn greatly suppresses 
loop- induced contributions to their masses |2^, as well 
as possible undesirable light-field-mediated effects. 



Let 1715^ = ?(^)- As that the field $ is a gauge singlet, 
we may allow it to travel in various dimensions. It is then 
reasonable to assume that its coupling to x depends on 
the volume of extra dimensions. We will adopt the ansatz 

with Vn the volume of the n extra dimensions where <i> 
can propagate. From our previous discussion, we may 
expect that Sm^^ ^ — ms^/lO, which yields SVn/Vn ^ 
0.1, ignoring factors of 0{n). Hence, we require that the 
volume of some extra dimensions change by about 10% 
over a Hubble time. 

We will not speculate on the nature of the poten- 
tial that gives rise to the aforementioned slow evolution. 
However, we point out that the initial change of ^ can- 
not be too fast, as it could yield unwanted consequences. 
First of all, if this change takes place too quickly over the 
time scales of freeze out process, one may end up with 
insufficient amounts of DM. To avoid this problem, we 
then require that at time tfo ^ Afp/Tj^, where T/o is 
the freeze out temperature, <C 0.1. We also require 
that the predictions of Big Bang Nucleosynthesis (BBN) 
be largely unperturbed. This implies that the kinetic en- 
ergy in the radion (5tr)^ must not be the dominant en- 
ergy density, so that a radiation dominated cosmological 
scenario persists (see Ref. for an alternative setup). 
Assuming a typical BBN temperature of T ~ MeV, dur- 
ing this era {t I s), we then roughly require 

9tr<MeV^ (BBN). (10) 

Here, we note that there is a map between the radion, 
corresponding to the volume of the extra dimensions, and 
a Brans-Dicke (BD) field; see for example Ref. [Q. The 
BD parameter ui 1/N — 1, where N is the total num- 
ber of extra dimensions and we have n < N . The agree- 
ment of various observations with standard cosmology, 
governed by general relativity, then places constraints 
on the value of uj. Solar system measurements require 
uj > 40000 today However, there are also stringent 
bounds that come from WMAP observations of the CMB 
and large scale structure, suggesting that |w| > 100 [p4| . 
Clearly, this is in conflict with having a dynamical radion 
field that varies over a significant fraction of Hubble time. 
To go further, we then consider a few possible scenarios. 

Since we are only interested in shifting the couplings 
of X to the SM after freeze out, it is, strictly speaking, 
unnecessary to assume any time variation long after this 
era. In generic WIMP scenarios, Tfo ~ 10 GeV, corre- 
sponding to t/o ~ 10^^ s. This is well before the CMB 
decoupling era {t ~ 10^ yr), and even that of BBN. How- 
ever, only signals related to extra dimensions, and not 
cosmological signatures of time variation, would likely be 
observable in this case. Even though this is a consistent 
option, we will not pursue it further here. 
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Another possibility is suggested upon noting that the 
BD constraints obtained in Ref.[|^ are in fact mostly a 
result of constraints on the time variation of Newton's 
constant Gn- This variation, in our setup, is a conse- 
quence of change in the volume of the extra dimensions. 
Note that we are not necessarily assuming that ^ depends 
on the total volume of the extra dimensions, which is gov- 
erned by the radion potential; we only assume that the 
coupling ^ depends on some of the radii. If some radii 
change in a way that ensures the volume of the extra di- 
mensions remains constant over cosmological times 
we may bypass the constraints on variation of Gn and 
hence the BD parameter uj. This will correspond to a 
Kasner-like behavior for the cosmology, where some di- 
rections expand while others contract. We will assume 
that such a solution can be constructed, leading to a con- 
stant extra-dimensional volume, while allowing some of 
the radii to vary with time. As long as ^ depends on 
a subset of such radii, we will achieve a situation where 
^, and therefore the localization of Xj can change over 
times of 0{Hq^), without an appreciable change in Gn- 
This could lead to models that are unencumbered by the 
bounds on uj, as assumed in the following. 

The basic requirements of our scenario leave a wide 
range of possibilities open to model building. As an alter- 
native example, we will outline a setup in which x starts 
out baryophilic, having C(l) couplings to the quark sec- 
tor. Over cosmological times, the evolution of ^ changes 
the localization of x, moving it towards leptons. At late 
times, X ends up being leptophilic and baryophobic. 

The above possibility, in fact, fits very well with one of 
the original motivations for fermion localization in extra 
dimensional models, that is the geometric suppression of 
proton decay p2[ , by separating leptons and quarks by a 
sufficient distance over the width of the wall. Even if the 
cutoff of the theory is Mp ^ 1 TcV, one can show that a 
separation s ~ 10//i will suffice to get a safely long-lived 
proton If X mainly couples to quarks initially, we 

need a new state with the right quantum numbers to pro- 
vide the necessary annihilation cross section. For exam- 
ple, a new scalar 9 with the conjugate quantum numbers 
of SU{2) singlet ur quarks will in principle be sufficient. 
As for (f) in the previous case, we need to assume that 9 
is charged under the same Z2 symmetry as x, (9) = 0, 
and mg > m^. In analogy with Eq. (||), 9 couples via 

02^ Xg^xOuR, (11) 

where Xg^fK is a 5D Yukawa coupling. 

Given a separation s ~ 10//i of leptons and quarks 
required to protect the proton from fast decays, we now 
require dm^^ ~ 20/i, after freeze out. Hence, some of the 
estimates from the above discussion will be modified here, 
but not enough to change the qualitative conclusions of 
the previous case. As the field x moves toward its cur- 
rent location, one can expect its couplings to various SM 



fields to rise and fall with time, in a model dependent 
manner. However, with separated quarks and leptons 
in well-motivated models, a few general features may be 
expected to arise. 

If the evolution of couplings takes place over cosmolog- 
ical times, we may start with a population of DM that 
mostly annihilates into quarks. This will entail a hard 
photon signal component coming from i:^ —> 77 at high 
rcdshifts. To get a rough estimate of the possible signal, 
let us assume that the photons originate from a galaxy 
cluster of size ~ 1 Mpc at a distance Dc from the earth. 
We take a DM number density ^ 10^^ cm^^ and an 
annihilation cross section (aav) ~ 10"^'^ cm'^/s, similar 
in size to the implied leptonic cross section at the present 
epoch. The annihilation rate into qq is then given by 

Ta l^nliaav) ^lO^U-\ (12) 

Given that at 0(100) GeV energies we can expect an 
0(10) multiphcity of neutral pions 0, the flux of pho- 
tons at earth is approximately given by 

/,^WZ^.^^10-(^^)%m-s-. (13) 

Hence, if the quark annihilation channel was dominant 
around 10^ years ago, Fermi gamma-ray space telescope 
may be able to detect |^ the photons from tt" 77 at 
the corresponding redshift. 

As X enters a 'fermionic desert' between the loci of 
quarks and leptons, its major couplings are governed by 
higher dimensional operators that are suppressed by pow- 
ers of the cutoff scale Mp ^ A"-'^ . Hence, one can expect 
a drop in the efficiency of DM annihilation over interme- 
diate red-shifts. Once x moves closer to the locus of 
leptons along the extra dimension, its couplings become 
leptophilic and larger, leading to enhanced DM annihila- 
tion, dominated by lepton final states, at lower redshifts. 

In this paper, we proposed that the observed cosmic 
ray excess in positrons, and not anti-protons, can be due 
to the annihilation of DM with time-varying couplings 
that become more leptophilic after freeze out. A possible 
consequence of this scenario is a red-shift dependent anni- 
hilation signal. We discussed how these features may be 
incorporated in a model where the SM fermions and the 
DM fermion x are localized along an extra dimension. 
If the localization of x moderately changes after freeze 
out, it could lead to enhanced coupling to electrons and 
hence explain the recent data. Depending on details, the 
annihilation at earlier epochs and high redshifts could in 
fact be dominated by quark final states, resulting in hard 
photons from 7r° decays. 

Note added: After this work was posted, a new data 
set from the Fermi collaboration on the e"*" -I- e~ cosmic 
ray spectrum p6[ | was released. The new measurements 
do not confirm the excess flux around 500 GeV reported 
by the ATIC collaboration g , but they still allow a DM 
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interpretation. In particular, a DM candidate with a 
mass of a few 100 GeV, leptophilic couplings, and an 
enhanced annihilation cross section at the present epoch 
remains a possibility pTj for which our work proposes a 
time- varying extra dimensional scenario. 

We thank Mark Trodden for discussions. This work 
was supported by the United States Department of En- 
ergy under Grant Contract DE-ACG2-98CH1G886. 
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